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Orientation  of  Nonspherical  Aerosol  Particles 


ftOL  ASETRACT  omppv 
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The  orientation  distribution  of  non-spherical  particles  within  a  cloud  was 
studied.  Theoretically,  the  subject  was  treated  for  a  general  laminar  flow 
of  intermediate  as  well  as  weak  and  strong  character,  and  for  a  turbulent 
flow  whose  Kolmogoroff  scale  was  much  larger  than  the  particles  themselves. 
The  laminar  flow  served  as  one  realization  for  the  turbulent  field  within 
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'  the  cloud. 

It  is  found  that  the  orientation  distribution  function  in  a  laminar  medium  is 
generally  structured  with  a  preferred  orientation  while  that  of  the  turbulent 
medium  is  random  on  the  average  but  with  a  significant  dispersion. 

Experimentally,  a  horizontal  wind  tunnel  was  erected  for  orientation  measure¬ 
ments;  it  was  equipped  with  a  dedicated  hot-wire  computerized  system  for  tur¬ 
bulence  characterization,  which  was  shown  to  conform  to  the  determinations  of 
Laufer.  Light  scattering  signals  were  used  as  the  orientation-gauge. 


I.  INTRODUCTION 


The  orientation  distribution  function  F  of  nonspherical  particles 
is  of  great  importance  for  the  evaluation  of  system  properties  with  a 
directional  dependence.  In  the  transport  mechanics  field  of  interest, 
the  mean  settling  velocity  of  a  cloud  of  orientable  particles  (1),  and 
the  particle-ensemble  orientational 1y  averaged  diffusion  tensor  (2-^' 
serve  as  examples  of  a  property  determined  by  F.  Also,  the  optical- 
orientation  interaction  in  a  system  of  morphologically  unisotropic 
particles  (5-7)  constitutes  another  case  of  this  behavior.  The  orien¬ 
tation  distribution  of  nonspherical  particles  or  molecules  in  solution, 
which  is  essentially  the  same  phenomenon,  has  been  a  major  subject  of 
interest  in  rheology  for  a  long  time. 

In  essence,  the  orientation  of  particles  in  space  depends  on  the  inter¬ 
play  between  a  randomizing  action  such  as  that  of  rotational  Brownian 
motion  and  microturbulence,  and  an  orienting  action  arising  from  fluid 
dynamic  field  gradients  and  external  torques j  likewise,  the  orientation  can 
be  affected  by  a  (physical)  concentration  source  due  to  translational 
diffusion  of  the  particles  (8). 


In  the  frequently  encountered  cases  characterized  by  the  absence  of  ex¬ 
ternal  orienting  forces  and  by  a  (concentration)  source  whose  contribution 
can  be  ignored  due  to  the  small  value  of  the  ratio  between  a  typical 
particle  size  and  system's  dimension  (9),  and  for  laminar  flows,  the  ori¬ 
entation  distribution  function  F  will  be  determined  only  by  the  flow 
gradient  and  particles'  rotational  Brownian  motion.  The  relative  weight 
of  these  factors  is  expressed  in  the  rotational  Peclet  number  «  defined  by 
(X  r W«/D„  where  Wo  is  a  typical  gradient  and  Do  is  a  typical  rotational 
diffusion  coefficient.  Essentially,  the  function  F  is  a  flow  property. 


Studies  on  the  nature  of  the  distribution  function  have  been  carried  out 
for  cases  of  spheroidal  particles  immersed  in  a  special  though  represen-  ^ 


□ 
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tative  in  certain  situations  flow  like  that  of  a  simple  shear  (10,11) 
and  an  irrotational  extensional  one  (12).  Also,  investigations  have  been 
performed  for  the  asymptotic  situations  of  a  very  strong  field,  <x  >_>  1 ,  and:zz 
a  very  weak  field  ot  «  1,  and  nearly  spherical  particles  (8,13).  _ 
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The  interest  in  rheology  was  mostly  in  the  second  moment  of  the 
particle  orientation  vector  p)  so,  for  the  intermediate  cases  of 
it  was  sufficient  (14)  to  use  a  closure  interpolation  technique  based 
on  the  solution  for  the  asymptotic  situations, and  construct  approxima¬ 
te  constitutive  equations  for  rather  than  find  the  function  F  it¬ 
self.  The  derivation  of  the  latter  function,  from  which  the  moments  f 
any  orientation-dependent  quantity  can  be  deduced,  is  very  desriable 
but  extremely  complex  for  intermediate  values  of  a  (14);  it  has  not  been 
achieved  up  to  now  for  these  m -cases  and  for  a  general  shear  flow. 

AIM  OF  STUDY 
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Since  many  aerosol  systems  in  closed  regions  involve  intermediate  Peclet 
numbers  and  multidimensional  flows  where  all  nine  components  of  the  gra¬ 
dient  tensor  U  may  be  present,  and  as  many  of  the  systems  contain  particles 
varying  in  their  typical  aspect  ratio-  over  a  wide  range,  it  was  worthwhile 
to  undertake  a  study  aimed  to  obtain  the  orientation  distribution  function 
F  for  such  situations.  In  general,  we  planned  to  conduct  the  research 
according  to  the  accepted  picture  of  a  turbulent  uni  (or  multi) -phase 
medium.  In  this  picture,  the  turbulent  field  average  properties  are  assu¬ 
med  to  be  created  by  the  contribution  of  (theoretically)  an  infinite 
number  of  the  so-called  "realizations"  in  each  of  which  the  motion  of  a  fluid 
element  or  a  particle  at  every  space-time  point  is  uniquely  defined. 


Thus,  it  was  natural  to  divide  our  study  into  the  following  three  parts: 

1.  A  (theoretical)  analysis  of  the  orientation  behavior  of  model 

nonspherical  aerosol  particles  within  one  "realization" 
which  is  characterized  by  a  laminar  medium,  in  essence. 


2.  A  (theoretical)  analysis  of  the  orientation  behavior  of  the 
above  particles  within  a  turbulent  field,  with  a  special  emphasis 
on  the  lower  atmospheric  boundary  layer  and  a  turbulent  jet, 

3.  An  experimental  test  of  the  theoretical  conclusions  by  some  effect 
such  as  light-scattering  . 

In  this  stage  of  the  study,  it  was  reasonable  to  concentrate  on  the  case 
of  axisymmetric  particles  such  as  spheroids  with  an  aspect  ratio  R(=a,/<3i) 
of  0  ■<  R.  oo  . 
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Needle-like  prolate  particles  may  simulate  asbestos  fibers  and 
straight-chain  aggregates  while  disk-like  oblate  particles  may 
represent  mineral  platelets,  thus,  though  the  problem  of  mode¬ 
lling  of  real  (geometrically  impure)  shapes  of  nonspherical  aero¬ 
sol  particles  by  regular  forms  is  still  unsolved  in  general,  the 
above  axisymmetric  shapes  can  serve  as  useful  approximations. 

In  passing,  it  should  be  stressed  that  the  herewith  reported 
research  was  devoted  to  such  a  case  in  which  the  particles  could 
be  considered  as  moving  in  a  continuous  fluid.  The  other  cases 
such  as  that  of  particles  comparable  in  (typical)  size  to  the  mean 
free-path  of  the  (gaseous)  medium  is  not  amneable  as  yet  to  a 
(sufficiently)  exact  analysis. 


II  .  THE  LAMINAR,  ONE  "REALIZATION"  CASE 
1.  Mathematical  Formulation  of  the  Problem 

Consider  a  system  of  small  spheroidal  particles  embedded  in  a  general 
shear  laminar  flow  u  which  can  be  given  in  their  vicinity  by  ^ 


u(x,t)  =  W(tMf. 


(2.1) 


The  distribution  function  F  is  described,  in  orientation  space  and  with 
the  justifiable  neglect  of  a  concentration  source  (see  above),  by  the 
known  Fokker-Planck  equation 


^  ♦  V.(Fw  -D  VF)  =  O 

W  V 

where  is  particle’s  rotational  velocity. 


(2.2) 


In  the  following  text  the  components  of  u>  for  rotational  ellipsoids 
in  the  general  laminar  flow  of  (2.1)  are  given  in  part  a  ,  which  can  be 
viewed  as  a  preliminary  treatment;  then,  the  solution  of  Eq.  (2.2)  for 
such  a  flow  is  presented  in  party,  and  numerical  experiments  for 
interesting  cases  are  brought  out  in  part  4  . 

Vor  nomenclature  see  APPENDIX  A  . 
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2.  The  Rotational  Velocity  Components  in  A  general  Shear  Flow 


In  order  to  solve  the  Fokker-Planck  equation  for  F  in  a  general  flow 
field,  one  has  to  know  the  dependence  of  the  rotational  velocity  CJ  on  W. 


The  components  U)  for  ellipsoidal  particles  in  a  creeping  flow 
were  calculated  by  Jeffery  (15)  in  general  but  were  specifically  given 
by  him  only  for  the  rotational  bodies  in  a  simple  shear  flow.  These  com¬ 
ponents  were  later  obtained  by  Hinch  and  Leal  (16)  for  any  ellipsoid 
though  also  only  for  the  above  kind  of  flow.  On  the  other  hand,  Takserman- 
Krozer  and  Ziabicki  (12)  presented  expressions  for  irrotational  flows  which 
hold  however  solely  for  the  ax i symmetric  case.  Thus  to  find  the  components 
of  cu  in  the  general  shear  flow  we  had  to  proceed  in  the  following  way: 

Let  us  take  a  rectangular  coordinate  system  X|^  fixed  in  space  and  a  (rectan¬ 
gular)  one  x'^  locked  in  the  particle  with  Euler  angles  d,  ^ such  that 
0  is  the  angle  between  x'and.  is  the  angle  between  the  planes  x, 
and  x,x,'  and  is  the  angle  between  the  planes  and  x'x^ 

(Fig.l).  Then  the  transformation  matrix  between  the  two  systems  x*  and  x^ 
will  be  given  by 


ce 

sOcjP 

sBs\p 

B  = 

■heecjpcf 

c^Sif 

4-cBifcy/r 

-cQcfsy/f 

(  2.3) 

Kinematically  the  components  of 

(jj  in  the 

system  are  related  to  the 

Euler  angles  and  their  temporal  change  by 

w,  =V^c0  . 

=  -  >jps9c\l^, 

Qc-d/  +  *fs6 

^  'f'  T  r  ^  ^  ^ 

where  c0  is  cose  ,  s9  is  sin0  and  similarly  for  iP  and  ^  j 

likewise,  the  components  for  sufficiently  small  ellipsoids  and  with  the 
neglect  of  torques  are  expressed  by  Jeffery's  general  relationship  (15) 


(2.5) 


where  i,j.  ^  — ►  1.2,3  are  cyclic  permutation  indices  ,a-,aj,a*  are 

the  semi-axes  of  the  ellipsoids  and  the  vorticities  Vj^-  and  the  strain 

coefficients  S--  are  defined  by 
•■j 

for  the  assumed  gradient  field  in  which 

K'i = i  (2-^) 

J 


(2.7) 


In  this  case  of  small  ellipsoids, Eq.  (2.5)  should  hold  even  for  non-uniform 
slow  motions  of  the  particles  if  the  typical  time  of  change  of  H(t)  of(2.1) 
is  large  enough  in  comparison  with  the  relaxation  time  of  the  latter. 


By  the  similarity  transformation  from  x,,  to  x^  we  have 


Hi  = 


(2.8) 


where  the  repeated  index  implies  summation  . 

,  /’'■> 

/ 

'<«  / 


Fig.l.  The  Euler  angles  and  two  systems 
of  coordinates. 


So  Eqs.  (2. 4), (2. 5)  and  (2.8)  for  rotational  ellipsoids  with  lead 

to  the  final  expressions. 

^  =  V,j/2  *  ctge(E„c)8-£„sW  +  >(^S„c^Si’-  Q,\1') . 
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®  '  G.,„  e>@ - S„ ,J  s’e*  I  S^e Q,c,»- *  I S„ s^ip).  ” 


Here,  for  brevity  sake,  we  denote 
^ik  ~  Z^^ik 

Q,=  i(w,,-W„)  ,  Qj  =  5(w„»w,,-2v><,), 

''i*  '  '^ik  ■  I  Sj,  =  ♦  W'^ ;  , 

W,=iiii  • 

‘k  JX,  ' 

c"*e=cos"’0,  s„0  =  sinme,  s'^e  =  sin"'e 


(2.10) 


and  similarly  for  j  >  is  defined  by  X  «  (k‘-1)/(R^*^1  )  . 

The  special  cases  of  the  previously  studied  flows  and  axlsynmetrlc  particles  are 
readily  deduced  from  the  above  equations.  Thus,  Jefferey's  frequently  used  solution 
(15)  for  a  Couette  flow  (0,0,qa^ ) follows  from  (2.7)-(2.8)  as 


a  -  > 


Likewise,  for  the  Irrotatlonal  flow  characterized  by 


(2.11) 


K.  =  '^55  “-'"'..A  =  -<f/i 


(2.12) 


(2.13) 


and  all  other  V/|^  equal  to  zero,  we  get 

>»-o,  e  *  J  QjS,e  = -|;Xqs,e, 


which  are  in  agreement  with  Takserman-Krozer  and  Zibabicki  (12). 
3.  The  Solution  for  the  Orientational  Distribution  Function 


Let  us  designate  by  a  velocity  gradient  typical  of  the  considered 

problem,  and  introduce  the  dimensionless  parameters 


t  =  wt . 
o 


(2.14) 


then,  the  diffusion  equation  (2.2)  together  with  (2.4) ,(2.9)and  (2.14)  can 
be  written  as 


iL  ^  (Fesine)  -r  i-  (Fjp)  =  AF/oc 

9t  Sine  30  ^  df  ' 

where  the  rotational  Peclet  number  o<  is  defined  as  before  by  ot  = 

and  0  now  is  the  mid-diameter  rotational  diffusion  coefficient;  the  sign  ~ 

is  henceforth  dropped  out  for  simplicity. 

The  normalization  condition  is  taken  as  usual  to  be 

^^^J^Fsinededf  =  1  .  (2,16) 

According  to  the  theory  of  Fourier  series,  one  can  obtain  the  solution  of  (2.15) 
in  the  form 

F  =  Z  f  (2.17) 

ft-o  m-o  ^  i  'WO/'-  on  »n'  in'  '  m' J  n 


as  long  as  the  equation  for  F  has  a  Laplacian  in  6  and  .  Such  a  function, 
which  is  twice  differentiable  with  regard  to  these  variables,  can  be  uniformly 
approximated  at  every  moment  t  by  a  finite  series  of  the  orthogonal  spherical 
harmonics. 


Inserting  (2.17)  in  (2.15)  ftnd  taking  •  into  account  the  syinnetry  properties 
of  F  based  on  physical  reasoning, one  gets  for 


dA^„  (t) 


n(n-t-1 ) 


**-♦■1  m*2  1 

-  r  I  £  (1 

j‘=n~i  psm-i  in 


■s™.) 


(2.18) 


where  k  =  0.1,  n  =  2,4  •  .  .  ,  and 


m  -  0,1  ,  2, 


pm  pm 

The  expressions  for  the  coefficients  and  appearing  in  (2.18)  are 

given  in  Table  1.  The  vorticities  ^  the  strains  S^pt  ^  and 

are  defined  by  (2.10);  the  indices  j  and  h  are  even.  The  numerical 
coefficients  dj^  can  be  calculated  from  Table  II  according  to  the  position 
indicated  in  the  last  column  of  Table  I.  In  Table  II  =  i/f2(2n*l)(2n  4l)J  , 

e^- l/[2(2n-i)(lri  +  l)]  ^  and  e ^  =  l/[2(2n+ 1  )(2n  +  i)]  j 

all  other  coefficients  including  those  of  negative  indices  and  those 

related  to  index  combinations  not  appearing  in  the  table, are  zero.  The  corres- 
ponding  Dij„  andC^  jn  coefficients  fall  out  too.  The  value  of  A^^-l/ln 
can  be  obtained  from  the  normalization  condition  (2.16)  and  =0  since 

mi^=  0  for  m  =  0.  Here  it  should  be  noted  that  for  the  steady-state  situation, 
when  dF/dt  =0  ,  the  r.h.s.  of  the  general  equations  (2.18)  provides 


recurrence  relations  for 
studied  special  cases. 


and^^^  which  include  those  of  the  previously 


TABLE  I  pim 

The  Coefficients  D,  .  and  C,  • 


1  m  ±  I  I.  2.  1- 

2  m±2  0.2.4. 


LocMNNi  of  dp*  im 
TaMc  II 


X4CI2<?.(*  -  *-«)  +  QMroJ  +  <#»  -  m)X<«rS,2/2  Ifp  =  m  -  0  III 

Else  II 


Now,  just  Compare  our  solution  with  some  of  the  former  studies,  let  us  take  the 
stationary  gradient  field  with  perpendicular  and  parallel  components 


r 


4S3 
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(2.19) 


which  includes  as  special  cases  simple  shear  and  extenslonal  flows.  Considering 
the  field  of  (2. 19), small  enough  >  values*  and  truncating  (2.17)  and  (2.18) 
at  Ki  =  2  .one  obtains  for  the  first  terms  of  the  steady-state  solution 

as  expanded  In  A 

P  -  ±  /l-  .  3As^e 

"4.w|  2f3A-i-AQ2)  2(36%PT^) 

Thus,  for  the  Couette  flow  {0,0, <JX^)  .when  only  Is  not  equal  to  zero, 

the  solution  described  by  (2.20)  coincides  with  that  of  Peterlln  (10)  In  the 
first  order  of  A  If  one  takes  Into  account  the  convention  difference  of  signs 
and  puts  In  (2.14)  =  .  Also  In  the  case  of  the  irrotatlonal  flow  of 

(2. 12), we  get  by  an  expansion  In  ^o<  an  accord  with  Takserman  Krozer 
and  ZlabKki  (12). 

These  agreements  could  be  anticipated  on  the  grounds  that  the  corresponding 
coefficients 


TABLE  II 

mm 

The  Numerical  Coefficients  ^rtn' 


n  -  2  m  -  I  -»i(«  -  2)  n  ("  +  "•  “  •)  m  -  2  ^  (w  -  2)  H  (n  +  m 

(-«  *  )•« 


-  i)  0 


n  m  -  I  ie^lm  -  l)/2 


2^1  * 

« -  2  -z"  n 

*  .-I 

X  (w  +  in  -  I  +  I) 


-'i  n  ~ ') 


#VI("  +  0 


If  +  2  m  -  I  -e^H  +  3)  n  (if  -  m  +  0  m  -  2  -  y  (»  +  3)  H  (»  -  hi  +  i)  0  17  (ii  i) 

1*1  ^1.1  i.i 

11-2  III  +  I  ri(ii  -  2Xii  +  Iff)  Iff  +  2  ^iCh  -  2V2 

If  fii+l  3e^2m  l)/2  m  2  2ej2 

11  +  2  m+l  rrfii  +  3K»i  -  m  +  I)  m  +  2  +  3V2 


.•.y* 


►  ,  •  ,  «  .  1 
C-I'Cn 


pMi  pm 

0.  >Ci  Tables  I  and  II  are  in  harmony  with  those  of  the  recurrence 

relations  of  both  the  latter  studies.  Obviously »  for  arbitrary  prameters  \o<.t  and 
one  has  to  take  a  sufficiently  large  N  which  requires  the  use  of  the  computer. 

4.  Numerical  Results 

The  investigated  particles  were  long  prolate  ellipsoids  of  aspect  ratio  R<  10^ 
and  thin  oblate  spheroids  of  aspect  ratio  0.01  ^  R  4  0.1  which  can, 

respectively,  simulate  long  cylinders  and  thin  platelets.  The  average  characteristic 
dimension  of  these  particles  was  taken  to  be  cm;  so  their 

rotational  mid-diameter  diffusion  coefficient,  as  calculted  from  the  classical  expre- 
ssions(see  Gans  (17),  was  O  ~  isec"^  .  The  typical  flow  velocities  of  interest 
were  assumed  to  be  1-10  m  sec'^  and  the  typical  velocity  gradients  1-10  sec'\ 
all  of  which  correspond  to  Peclet  numbers  of  ot  1-10. 

We  considered  various  models  of  flow  which  enabled  us  to  stydy  possible  effects 
on  the  orientation  of  the  aerosol  particles  in  closed  environment  and  atmospheric 
conditions.  The  initial  orientation  of  the  particles  was  taken  as  random;  the  field 
could  be  either  stationary  or  nonstationary.  In  the  last  case,  the  components 
were  supposed  with  no  loss  of  generality  to  vary  with  time,  with  a  frequency  A 
and  initial  phase  >9  as  cos  (At 

The  system  (2.18)  was  treated  by  a  high-accuracy  differential  equations- 

solver  based  on  the  extrapolation  method  of  Burlish  and  Stoer  (18).  All  the  summa¬ 
tion  operations  in  the  r.h.s.  of  (2.18)  as  well  as  in  (2.17)  were  carried  out  with 
double  precision.  The  program  was  checked  against  the  solution  for  F  obtained  by 
Peterlin  (10)  and  Takserman-Krozer  and  Ziabicki  (12)  for  the  hbove  cited  special 
cases.  For  stationary  flows  and  sufficiently  large  t  when  a  steady  state  is  being 
reached,  our  solution  coincided  with  that  achieved  by  us  also  through  the  Gauss  eli¬ 
mination  method  applied  with  double  precision.  This  has  been  carried  out  to  test  the 
computational  compatibility  of  different  methods  in  treating  the  system  of  (2.18). 

A  fast  enough  convergence  of  the  solution  for  F  with  a  reasonable  accuracy  has  been 
achieved  as  a  rule.  The' calculations  were  carried  out  for  oc.  with  N  *4,6 
(n,N-even)  which  coincided  to  up  to  four  decimals  and  better  with  solutions  where 
a  higher  value  of  N  was  adopted;for  oc  *  10,  we  used  N  *  10-14  and  got  an  accuracy 
for  maximum  and  mean  values  of  better  than  1%  which  was  considered  quite  satisfactory 
for  our  purposes.  In  the  case  of  |3k|  ^  0.923  (0.1^  5)  or  some  fields 

of  smaller  values  of  N  could  be  taken. 


Of  the  (physical)  situations  which  were  considered,  we  bring  here  these 
of: 

a.  A  point  source 
and 

b.  A  circular  laminar  jet 


A  point  source  flow 

We  consider  here  a  flow  created  by  a  point  source,  where  the  velocity  is  given 
by  and  its  gradients  are 


W: 


(3) 


ik 


(2.21) 


with 


is  nondi mens ionali zed  by  a  typical 
length  is  the  radius  vector  and  a  constant.  is  conside¬ 

red  to  be  constant  in  the  immediate  neightboHiood  of  the  (small)  particles. 


The  orientation  distribution  function  F  for  such  a  flow  is  presented  in  Fig.  2 
as  F(9^f)  vs  6  for  y  =  o^tc/4.  j  oc  =  i;  R=50(X5:  1  )  5  and  R.  =  o.02(> 

,  a  platelet),  and  for  the  (arbitrary)  space  point  f 4.. 

At  this  point  we  have 

-.84.0  -.599  -.599 

-.599  .^20 

-.599  -.2fO  .f20  /  .  (2.22) 


7rom  Fig.  2  it  can  be  seen  that  ,  as  soon  as  the  components  of  decrease 

with  distance  and  become  mostly  negative,  the  location  of  the  maxima  of  F  for  cylin¬ 
drical  and  platelet-like  particles  change  to  values  of  6  and  ^  opposite  to  those 
of  the  former  cases  of  fields  and  .  In  Fig.  3  are  brought  out  curves 

of  the  maximal  values  of  F,  F_  ,  as  a  function  of  timet  for  various  particle 
aspect  ratios  0.02  <  R  ^  SO,  1  <  oc  ^  IO  ^  and  for  the  same  space 

points  as  of  Fig.  2.  The  positions  of  the  first  maxima  in  F  are  shown  too. 

An  important  point  to  be  noted  in  this  flow  case  is  that  the  orienting  action 
of  the  velocity  gradient  (2.22)  is  noticeably  stronger  for  platelets  (curves  1,2) 
than  for  cylinders  (curves  7,8)  at  oc  -  10  but  that  this  difference  is  weakened 
with  the  decrease  of  oc  . 


Fig.  2  .  The  orientation  distribution  function  F  vs0for  the 
point  source  gradient  of  (1.21)  ,  a  =1,R  -O.OZ 

(platelets)  and  R  =  50  (long  cylinders)  \  >^1\  X^  =  1,Xj  =  x^  =  O.V- 

A-f  =0;  - ,  =50:  i.  ya-o^2.  y*sn/v; - ,  R  =o.oa  :  3.  )P*o^ 
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1.  ot 

2.  oi 
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3.  ot  =  S’,  R  a  0.02 
4..  ot  ®  5,  R  aO.l 
S.  tK  =  1 ,  R  =0.0  2 
b.  oc  *  1 ,  R  » O.f 
01  s  fO,  R  »50 
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f  O.  ec  =  5,  R  =  1  O 
11.  «».  =  1,  R  «  50 
f2.  oc  a  1,  R  a  10  , 


(t  >  1; 

e^  =  3»rA 

Pm  -  «/V 

9m*  2-^/^ 

(t^t) 


Fig.  3 

with 


vt.  t  for  a  point  source  gradient  (2.2 

X3  =  0.4  . 


Fig.  4  A  perspective  view  of  particle  orientation  in  a  point  source  field 
(scheinatical)  . 


A  circular  laminar  .let 

The  velocities  within  a  circular  laminar  jet  are  given  by  Schlichting  (19)  as 


1 


(2.23) 


where  ^  /(f.)  .Here  V  is  the  kinematic 

viscosity, -y  is  determined  by  the  jet  momentum  J:  ~n/>g*vj ^  density ,  and  r  is  the 
radial  distance;  t-  and  are  nondimensionalized  later  on  by  again. 


For  this  flow,  the  components  of  Wf*^,r)  are 
= 


(2.U) 


with  ^=2-jr^y/[xj(u^/^fy,they  are  taken  to  be  constant  in  the  vicinity 
of  the  (small)  particles  as  in  the  point  source  case. 


The  dependence  of  F  on  6  at  'P  -0, 11/4 
for  the  indicative  case  of  <>(=1,7=1 


is  presented  in  Fig. 5 
various  particle  aspect  ratios  R, 


the  space  point  »1 

%  =0.5,  X  =  X  , 

and  = 

Here  the  velocity  gradients  which  decrease  with 

the  distance  too,  are  (in 

Cartesian  coordinates  ) 

/-.678  -.295 

-.295  \ 

I  -.239  .359 

-.0?7  ) 

\-.259  -.0^7 

.359  / 

(2.25) 

and  thus  one  can  conjecture  that  the  character  of  the  orientation  distribution 
function  of  platelets  and  cylinders  would  be  similar  to  that  of  the  point 

source  case. 


A  perspective,  schematical  illustration  of  the  orientation  in  the  jet's  field 
is  shown  in  Fig.  6  which  is  constructed  with  the  consideration  on  the  jet's 
axial  symmetry. 


Fig.  5.  The  orientation  distribution  function  F  vs.  d  in  a  flow  of 
axi symmetric  laminar  jet  (z.24)-(£.2.f)  at  X,  =1^  %  =0.5^ 

for  platelets  and  cylinders  at  R  =  0.01,  0.2,5,100  ;  ^ t>2 

1.  R*O.Of, /»=Oi  3.  R=0.2,/»  =  0;  5.  R  =  /oo,  y>  =  o;  l-.R-S,f  =  Oi 

2. R-o.oi,f=rc/^  ,  A-  i  6.  =  i  «.  R=5, 


Fig.  6.  A  perspective  view  of  particle  orientation  in  a  laminar  jet  field 
(schematical) 

The  computer  code  for  the  (steady) laminar,  one"i<ealization''  case  is  deposited 
with  the  CDRC,  Aberdeen  Proving  Ground,  HD.,  c/o  Dr.  J.  Buttiger. 


5.  Interim  Conclusions 


Obviously,  we  did  not  intend  in. this  study  to  cover  every  significant  case  of 
a  particle  assembly;  however,  the  type  of  particles  and  flow  fields  covered  by  us 
and  given  here  is  very  conductive  to  the  prediction  of  the  behavior  of  real  and 
important  systems.  The  following  conclusions  can  be  drawn  out: 

(i)  The  method  put  forth  in  this  study  provides  a  useful  and  convenient  tool  in 
the  structure  analysis  of  the  orientation  distribution  function  F  for  a  wide  class 

of  particle  parameters  and  flow  configurations.  In  practical  use,  the  obtained  solution 
converged  fast  enough. 

(ii)  The  solution  for  F  achieved  here  includes  as  special  cases  those  of  previous 
studies. 

(iii)  The  character  of  F  has  some  regularities  distinguished  by  perceptible  maxima 
whose  values  and  locations  can  be  found  in  a  general  shear  flow;  these  maxima  will 
give  clear  indication  of  a  preferred  orientation. 

(iv)  The  time  evolution  of  F,  the  time  relaxation,  and,  correspondingly,  other 
integral  characteristics  of  it  can  be  determined  now  with  sufficient  accuracy. 

(v)  The  positions  of  the  maxima  of  F  are  established  fairly  fast  and  do  not  change 
appreciably  with  time. 

(vi)  Significant  differences  exist  in  the  orientational  effect  regarding  the 
maximum  values  of  F  for  platelets  and  long  cylinders  in  different  flow  situations; 
the  effect  is  stronger  for  platelets  than  for  cylinders  in  velocity  gradient  fields 
which  decrease  with  distance. 


APPENDIX  A  ;  NOMENCLATURE 


semi-axes  of  ellipsoidal  particle  (a,  .lies  along  x,  ) 

A  transformation  matrix  with  component 

coefficients  in  expansion  of  F,  Eq.  (2.»’») 

B  transformation  matrix  with  components  Eq.  (2.3);  .bT,^  components 

of  inverse  matrix 

c^O  cos  mO  .and  the  same  for  and  ;  C*'®  cos^0  .and  the 

.  am  pm 

same  for  ^  and  coefficients  in  Eqs.  (2.i8) 

and  Table  1  ' 

0  rotational  diffusion  tensor  of  the  particle;  D  rotational  diffusion 

coefficient  around  a  mid-diameter 

defined  in  text,  Eq.  (l.itf) 

®i*®i**a  coefficients,  Table  II 

F  orientation  distribution  function  •  .  maximal  value  of  F 

^  uikt  defined  in  text.  Eq.  (2*I0; 

C 

running  indices 

in,  m' 

n,n',  p 

J  jet's  momentum 

N  upper  limit  of  n  in  Eqs.  {2.f*»)and  (2.  U) 

m 

associated  Legendre  polynomial 
q  velocity  gradient 

constant 

functions  of  fluid  flow.  Eq.  (  2.io) 

R  particle's  aspect  ratio  R  =  a^/a^ 

r-g  typical  length 

strain  coefficients  referred  to  coordinate  system  coefficients 


referred  to  system  xj^ 


6  sin  rw0  and  the  same  for  'P  and  ^ 

"'S  sink's  and  the  same  for  p  and 

t  time 

I  fluid  velocity;  components  in  coordinate  system  locked  into  the 
particle 

/ 

•j  vorticity  coefficients  referred  to  system  coefficients  referred 

to  system 

'  velocity  gradient  tensor  ;  components  of  tensor  referred  to 

system  components  referred  to  ,  typical  velocity 

gradient 

«i»  (R-1,2,3)  Cartesian  coordinates  fixed  in  space 

x^  coordinate  fixed  in  the  particle 

Greek  Letters 

cx  rotational  Peclet  number 

/9  initial  phase 

y  constant  determined  by  the  jet's  momentum  J 

S  Kronecker’s  delta 

Euler's  orientational  angles  ;  0^  , 
are  positions  of 

dimensionless  particle  parameter  defined  by  ^  0 

rotational  velocity  of  the  particles  with  components 
frequency 
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III  THE  TURBULENT  CLOUD  CASE 


1.  The  (Theoretical)  Models 


Having  performed  a  study  on  the  orientation  density  function  (o.d.f.)  of 
small  spheroidal  aerosol  particles  in  a  general  ^ield  of  an  arbitrary  strenrth, 
(1,2)  it  became  possible  to  treat  the  orientation  problem  in  a  turbulent  medium. 

To  this  end,  two  interconnected  physical  models  were  applied.  In  the  first. 

"The  Realizations  Model,"  it  was  assumed  that  the  turbulent  particle  field  con¬ 
stituted  an  ensemble  of  an  infinite  number  of  realizations,  j,  each  one  of  which 
is  characterized  by  one  set  of  the  o.d.f.,  F^'^^(x,t)  values.  The  latter  was  taken 
in  that  model  to  essentially  coincide  with  the  previously  found  solution  of  the 
Fokker-Planck  equation  (i^2)  in  the  field  of  the  realization 


|L  -  ^ 


0) 


,  (j> 

U)(  , 


R)  F  -  \7F  J  = 


(3.1) 


where  cuis  the  rotational  velocity  of  the  particles  and  ot  is  their  rota- 

/F*  ^ 

V/.  being  a  typical  component  of 

the  (fluid)  gradient  tensor,  an  effective  rotational  diffusion  coefficient 
and  R  the  particle  aspect  ratio  j  W.j^  is  a  gradient  component 
a  component  of  the  fluid  velocity).  The  use  of  the  solution  is  based  on  the 
estimate  that,  even  for  the  highest  (Kolmogoroff)  frequency  component  of  tij  , 
the  rotational  Reynolds  nuiMber  for  the  studied  particles  is  small  enough  to 
render  their  motion  (quasi) stationary. 


In  the  second, "The  Micro-turbulence  Model,"  a  relation  between  the  so-called 

IT 

turbulent  diffusion  coefficient  0^  of  the  particles  and  the  physical  characte¬ 
ristics  of  the  fluid  field  was  supplied.  This  coefficient  and  the  Brownian  ro¬ 
tational  diffusion  one,  0^  .compose  the  effective  diffusion  coefficent 

\  =  \  \  .  (3.2) 

Expressing  the  stochastic  quantities  of  the  system  as  F  =  F *F\j.  ^  u. 
and  u=u*u',  it  could  be  shown  that  the  realizations  average  of  F,F{ 

-  ill  ^ 

for  a  space-time  point  is  actually  F(V^jj,  0^,R) 

for  that  point  and  that 


^  i 


The  particles  considered  were  taken  to  be  much  smaller  than  the  Kolmogoroff  scale;  so, 
the  turbulent  rotational  diffusion  coefficient  itself,  ^  assumed,  in  the 
second  model,  to  depend  on  the  randomizing  action  of  the  turbulent  pressure  fluctu¬ 
ations  at  the  particles'  surface  which  arise  from  the  (Kolmogoroff)  micro- turbulence. 

From  dimensional  analysis  it  was  obtained  that 

1“  it 

a  (£/v)  (3.4) 

where  £  ,  the  turbulent  dissipation  energy,  is  given  for  a  homogenous  field 

by  (3) 

£  =  (3.5) 

Ik 


2.  Applications 

As  the  turbulent  gradient  tensor  is  not  completely  known  in  many  fluid 

systems,  it  was  necessary  in  the  present  study  either  to  supplement  the  missing 
data  by  the  results  of  numerical  simulations  or  (reasonable)guesses  or  both.  In 
cases  of  interest,  there  was  conducted  a  parametric  investigation  in  which  the 
effect  of  structure  changes  in  the  gradient  tensor  on  the  o.d.f.  was  tested. 

The  numerical  simulations  of  the  turbulent  field  were  carried  out  according  to 
the  method  of  Wang  and  Frost  (4)  where  u\  =u\(t)  is  found.  In  this 

method,  however,  only  the  three  components  V/'f  (i=  1,2.3)  could  be  acquired  due  to 
the  applicability  of  the  relationship  known  as  the  Taylor  hypothesis,  viz. 


St  ■  ‘  jx, 


(3.6) 


to  those  components  alone. 


Resul ts 


Tmo  physical  situations  were  studied: 

1.  The  near-ground  atmospheric  boundary  layer,  in  which  only 


~  Ln(1  +  Y  ) 


and 


'^3  Sz 


(3.10) 


(3.11) 


essentially  exist,  z  being  the  elevation  above  ground,  the  mean  horizontal 

velocity,  u*  the  friction  velocity,  k^von  Karman  constant  and  z^the  roughness 
height. 

2.  The  turbulent  round  free  jet,  for  which  the  Schlichting's  solution  for  the 
average  velocity  components  (7)  was  used.  In  this  solution,  the  (molecular)  kinema¬ 
tic  viscosity  v  is  replaced  by  a  virtual  (turbulent) kinematic  viscosity  6^ 
given  in  the  equation 

/  iA 

£^s0.t6lj  (3.12) 

where  J'  J'  is  the  kinematic  momentum  of  the  jet, 

is  its  half-width  and  is  the  average  fluid  velocity  along  its  axis 

(  C  X,  in  which  c  is  0.63  to  0.79,  as  experimentally  found). 

In  the  boundary  layer  situation,  both  the  case  of  a  "weak  turbulence",  where 
W  and  a  "strong  turbulence",  where  W.  »  W,  were 

00^  Q  O 

investigated. 

For  the  first  case,  numerical  simulations  were  employed  while  for  the  second 
one  Klebanoff's  experimental  data  (5)  supplemented  by  simulations  was  applied. 

Also,  in  the  latter  case,  the  version 


was  taken.  The  calculations  in  both  cases  were  carried  out  with  ♦=0.15  cm^ec.,  S 
C  thickness  of  the  atmospheric  boundary  layer)  =  I0*m^  Zp  =0.1 -m,  k^=0.4 

Dg  =  Isec'*  and  e  deduced  from  Eq.  (5).  The  calculation  time  for  F  was 
greater  than  the  relaxation  time  of  the  particles,  and  A9=  =  n/ii  or  'ft/xif. 

Ja 

In  the  turbulent  jet  situation,  the  ratio  between  the  various  values  of 
was  obtained  (6)  through 

72  72 

C^W..  =  c=lc)  (3.14) 

where  c^- 1  *  exp(.2oor^)^  r  being  the  radial  distance  within  the  jet. 

For  both  situations,  the  absolute  values  of  were  deduced  by 

normalization  according  to  Eq.tJ.5).  Values  of  F,  and  Fm  ,  together  with  the 
deviations  of  F  for  fields  (2)4(3)  (above),  were  calculated  as  a  function  of  height 
above  ground  (and  hence  6  )  in  the  boundary  layer  situation,  the  aspect  ratio 
of  the  particles,  R  (R>i  for  fibers  and  R<  1  for  platelets),  and  the 
parameterization  factor  .  However,  only  the  following  typical  figures 

were  presented  here:  0.I2 

F 

O.K) 

aoe 

0.06 

0.04  _ _ 

Tr/4  ir/2  iwM  W 

_  9 

Fig.l.  Theaverage  o.d.f.,  F  ,vs.  ©  for  the 
atmospheric  boundary  layer,  "weak  turbu¬ 
lence"  ,  with  ; 

z  =  2.5  to  20m;(l )2.5m,(2)5m,(3)10m,(4)20m 

i,w;jf=Wo)  according  to  Eq.C3.ll)  ,  wj  =AV/,j. 

Solid  lines  are  simulation  results  (254  realizations,  usually);  points  are  values 
calculated  from  the  W. .  field. 


Fig. 2.  The  maximal  o.d.f.  ,  vs.  height  for  the  atmospheric  boundary 
layer,  "strong  turbulence",  R  =  SO;  e  decreases  with  height  accor> 
ding  to  Ball  (8).  (1)  w/jJ’  ,  (2)  (3)  K*”; 

lines  are  cases  where  terms  of  Eq.(3.7)  are  taken  with  opposite  signs 
points  relate  to  combination  of  Eq.  3.13b  (Line  1  essentially 
coincides  with  the  random  distribution  one). 


Fig.  3.  The  maximal  o.d.f.,  Fvn  vs.  particle  aspect  ratio  R  for  the 
atmospheric  boundary  layer ^"strong  turbulence". 

2  -3 

z»  5m,  6  ■  230  cm  /sec 

(1)  w[V  .(2) 

(Line  1  essentially  coincides  with  the  random  distribution  one). 


Fig.  4.  The  maximal  o.d.f.  F^miVS.  the  parameter  oc^  for  the 
atmospheric  boundary  layer;  "strong  turbulence". 

R  -50  - R*o.oa - •  z>  -2.5m;  £  =  5cm^/sec^. 

(1)  ,  (2)  W/J’  ,  <3)  V/jJ” 


K)-*  10*^  KT*  B  *0*  10*  10* 


Fig.  5.  The  maximal  o.d.f.,  vs.  particle  aspect 

ratio  R  for  a  free  turbulent  jet. 

Locations*,  x,  =2,  x^=  Xj  =  0.05  (see  ref.l), 

=  5  sec"^ ,  6  “  250  cm^/sec^ 

(i>,(2'U3V terms  of  Eqi3.7)  are  taken  with  opposite  si5iis. 
(Line  1  essentially  coincides  with  the  random  distribution  one) 


1.  The  average  o.d.f.  in  a  weak  turbulent  field  shows  structured  (preferred) 
orientation; 

2.  The  maximal  (and  average)  o.d.f,  in  a  strong  commonly  occurring,  turbulent 

field  of  the  average  realization  essentially  coincideiwith  the  random 

distribution  while  the  spread  of  the  values  of  the  function  between  fields 

quite  significant.  This  spread  of  values  may  have 

practical  connotations. 


3.  The  values  of  and  its  deviations  increase  with  the  parameter  ot  as 
expected. 

Free  turbulent  jet; 

Conclusions  2  and  3  of  the  former  situation  apply  here  too. 
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IV.  EXPERIMENTAL 


1 .Method 

The  experimental  method  consisted  of  a  determination  of  a  physical  effect 
from  which  the  orientation  density  functions  (o.d.f.)  F  of  the  particlescan 
be  deduced. 

Light  scattering  was  chosen  to  be  that  effect.  Thus,  the  method  adopted  by  us 
was  essentially  an  inversion  of  the  signal -integral 
iff  ITT  It 

7  =  /  I  f  // 

o  o  o  L  y 

drdL  sin  6d6  d<P 

(4.1) 

where  I  is  an  average  (measured)  light  intensity  at  an  angle  from 

the  initial  direction  of  its  propagation  ,  f(r)  and  g(0  are 

distribution  functions  of  characterizing  geometric  parameters  of  the  particles 
(r-radius,  1-length  for  cylinders),  and  ®  /  V'  Eulers' 

angles  (  >f/  disregarded  for  centrally  axi-symmetric  particles). 

Here,  as  one  did  expect  that  in  a  turbulent  field  particle  orientation  F  will _ 

fluctuate,  we  planned  to  also  determine  the  (central)  variance  of  tKe signal, 
given  by 


_  zn  xn  n 

”  O  o  O  I  y 

drdL  sinQdLQd'Pd 

(4.2) 

It  should  be  remarked  that  the  dependability  of  this  method  In  extracting 
a  sufficiently  accurate  values  of  the  o.d.f.  F  from  the  signal  integrals 
(4.1)  and  (4.2)  hinges  on  the  pre-knowledge  of  f  (r)  and  g(l).  So,  envisaging 

some  experimental  difficulties  in  this  direction*  we  launched  an  Independent 
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(and  financially  unsupported)  study  of  a  direct  determination  of  F. 

2.  Apparatus 
General  considerations 

The  experiments  were  conducted  in  the  fully-developed  turbulent  region 
of  a  straight,  horizontal  round  pipe  where  quiet  a  bulk  of  characterizing 
information  exists  (1).  Thus,  for  the  (round)  pipe  geometry^ one  has  that 
in  a  fully  developed  flow  (1) 

^  (4.3) 

/ 


Ur  .  =0, 


where  x  is  the  coordinate  along  the  average  flow,  and  U.^,U.r,U^ 
are  respectively  the  velocity  components  in  the  x,r  (the  radial)  and 
(the  azimuthal) directions. 

Also,  one  has 


u'rU'^=0  aU  T-  (4.4) 

where  the  prime  symbol  signifies  the  fluctmting  part  of  the  components. 
Specific 

A  schematic  diagram  of  the  apparatus  is  shown  in  Fig.  1- 


out 


Fig.l.  The  apparatus  (schematical ) 


in  which 


A  is  a  5  m  long,  5  cm.  i.d.  cylindrical  tunnel, 

B  is  an  aerosol  observation  section, 

C  is  Spurney's  fiber  aerosol  generator, 

D  is  a  light  source, 

E  is  (scattered)  light  detector(EMI  photomultiplier  type  9558)  and  read-out, 
^2  are  aerosol  sampling  ports. 

S  is  a  3  D  (or  2  0  X  -wire)sensor  (TSI,  Model  1295  -T1.5)  including: 

1.  CTA  supplier  and  linearizer  (0isa,56N21  with  56  06  Bridges) 

2.  Mini  computer  (Digital,  LSI/11/23  -*^,250  K  samples/sec) 

3,4,5  -  Diskette  drivefsw.^wrofa2s)»printer(FAClT)  monitor  (Visual  102). 

3.  Aerosol 

The  aerosol  used  consisted  of  glass  fibers  produced  by  the  Spumy 
generator  (Fig. 2a)  (2);  a  typical  size  distribution  of  the  fibers  is 
presented  in  Fig.  2b  (3) 

A 


D  ' - ’ 

Fig.  2a.  The  Spumy  generator  (schematical) 

F-container  of  fibers'  mass;  SE,RC=electromechanical  device 


Fig.  2b.  A  typical  size  distribution  of  the  glass  fibers,. 

S  4 

A  -  diameter  X  10  cm,  B  -  length  X  10  cm 


4.  Experimental  Procedure 

Setting  the  average  air  floif  through  the  tunnel  at  a  certain  value  befitting 
say,  an  average  Reynolds  number  of  10^>  the  turbulent  velocity  components  along 
the  X  and  r  coordinates  are  sampled  at  a  prescribed  rate  by  the  (computerized) 
system  1-4  (Flg.l)and  the  data  processed  for  .The 

necessary  gradients  (see  Part  III  of  this  report),  say 
can  be  deduced  from: 


1.  The  Taylor  hypothesis 


2.  The  incompressibility  demand, 

3.  Supplementary  values  taken  from  Laufer  (1)  for  the  dissipation  terms, 

4.  The  isotropic  conditions  assumed  to  hold  far  from  the  walls  of  the  flow 
tunnel  ( 1 ) . 

Then,  the  aerosol  generator  is  operated  and  aerosol  samples  taken  through 
portholes  Pp  P2  (FIG.1)  with  the  aid  of  a  (Cassela)  thermal  precipitator 
for  a  scanning  electron  microscope  analysis.  The  (aerosol)  collecting  surface 
within  the  precipitor  Is  wetted  before  sampling  to  obviate  the  need  for 
photograirmetry  (4). 

Finally,  scattered  light  signals  are  taken. 

Results  are  being  accumulated  now. 

V.  SUMMARY 

A  theoretical  method  has  been  developed  for  the  calculations  of  the 
orientation  density  function  of  elongated  and  flat  aerosol  particles  subjected 
to  the  most  general  laminar  flow  or  to  a  turbulent  field. 

Some  typical  situations  such  as  a  point  source  flow,  a  laminar  round  jet 
flow,  the  (lower) atmospheric  (turbulent)  boundary  layer  and  the  turbulent 
round  jet,  were  specially  treated. 

An  experimental  apparatus  essentially  consisting  of  an  horizontal  (pipe) 
wind  tunnel  has  been  constructed  to  check  the  theoretical  results  through 
the  effect  of  particle  orientation  on  light  scattering.  The  apparatus  , 
equipped  with  a  (dedicated)  very  fast  computerized  system,  is  functioning 
smoothly.  Scattered  light  readings  are  being  taken. 
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